Quantum synchronization and entanglement of two qubits 
coupled to a driven dissipative resonator 
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Using method of quantum trajectories we study the behavior of two identical or different supercon- 
ducting qubits coupled to a quantum dissipative driven resonator. Above a critical coupUng strength 
the qubit rotations become synchronized with the driving field phase and their evolution becomes 
entangled even if two qubits may significantly differ from one another. Such entangled qubits can 
radiate entangled photons that opens new opportunities for entangled wireless communication in a 
microwave range. 
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Recently, an impressive experimental progress has 
been reached in realization of a strong coupling regime 
of superconducting qubits with a microwave resonator 
[1, 2, 3, 4, 5]. The spectroscopy of exchange of one to 
few microwave photons with one [1, 3], two [2, 4] and even 
three [5] superconducting qubits has been demonstrated 
to be in agreement with the theoretical predictions of the 
Jaynes-Cummings and Tavis-Cummings models [6, 7, 8] 
even if certain nonlinear corrections have been visible. 
Thus the ideas of atomic physics and quantum optics 
find their promising implementations with superconduct- 
ing macroscopic circuits leading also to achievement of 
single artificial- atom lasing in a microwave range [9]. 

In comparison to quantum optics models [6, 7, 8] a new 
interesting element of superconducting qubits is the dissi- 
pative nature of coupled resonator which opens new per- 
spectives for quantum measurements [10]. At the same 
time in the strongly coupled regime the driven dissipa- 
tive resonator becomes effectively nonlinear due to inter- 
action with a qubit that leads to a number of interesting 
properties [11, 12, 13]. Among them is synchronization 
of qubit with a driven resonator phase [13] correspond- 
ing to a single artificial- atom lasing realized experimen- 
tally in [9] . The synchronization phenomenon has broad 
manifestations and applications in physics, engineering, 
social life and other sciences [14] but in the above case we 
have a striking example of quantum synchronization of a 
purely quantum qubit with a driven resonator having a 
semiclassical number of a few tens of photons [13]. This 
interesting phenomenon appears above a certain critical 
coupling threshold [13] and its investigation is now in 
progress [15]. 

Due to the experimental progress with two and three 
qubits [2, 4, 5] it is especially interesting to study the 
case of two qubits where an interplay of quantum syn- 
chronization and entanglement opens a new field of in- 
teresting questions. The properties of entanglement for 
two atoms (qubits) coupled to photons in a resonator has 
been studied recently in the frame of Tavis-Cummings 



model within the rotating wave approximation (RWA) 
[16, 17]. Compared to them we are interested in the case 
of dissipative driven resonator strongly coupled to qubits 
where RWA is not necessarily valid and where the effects 
of quantum synchronization between qubits and the res- 
onator is of primary importance. In addition we consider 
also the case of different qubits which is rather unusual 
for atoms but is very natural for superconducting qubits. 

In absence of dissipation the whole system is described 
by the Hamiltonian 

H = huoifi + 1/2) + nVticr^p /2 + nQ.2crf^ /2 (1) 
^ghujQ^cr^^ + cr^^){a + a^) + / cosut • (a + a^) 

where the first three terms represent photons in a res- 
onator and two qubits, _g-term gives the coupling between 
qubits and photons and the last term is the driving of 
resonator. In presence of dissipation the resonator dis- 
sipation rate is A and its quality factor is assumed to 
be Q = cJo/A ~ 100. The decay rate of qubits is sup- 
posed to be zero corresponding to a reachable experi- 
mental situation [1, 2, 3, 4, 5] where their decay rate 
is much smaller than A. The driving force amplitude is 
expressed as / = fiX^Jn^ where Up is a number of pho- 
tons in a resonator at the resonance uj = ujq when ^ = 0. 
The whole dissipative system is described by the master 
equation for the density matrix p which has the standard 
form [18]: 



p= - -[H,p] -^\{apd: -d:ap/2- pb:a/2). (2) 

The numerical simulations are done by the method of 
quantum trajectories [19] with the numerical parameters 
and techniques described in [13]. 

To analyze the system properties we determine the 
spectral density of driven qubits defined d^s S{v) = 

2 

/ (it exp{— zi/t}Tr(p(cri"^^ + cri^^)/2) . Its dependence 

on system parameters for identical and different qubits 
is shown in Fig. 1. At small couplings g the spectrum of 
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FIG. 1: (Color online) Spectral density S{h') of two driven 
qubits as a function of system parameters for identical 
Qi/cjo = 0.2/ ojQ — 1.2 (top row) and different Qi/cJo = 1-1, 
^2/(^0 = 1.2 (bottom row) qubits: (a,b) Up — 15, cj/cjq = 1; 
(c,d) g = 0.04, uj/ujQ = 1; (e,f) = 15, ^ = 0.04. Here and 
below A/cJo = 0.02. 



qubits S{v) shows the lines at the internal qubit frequen- 
cies VLi^2 but above a certain critical coupling strength 
g > Qc the quantum synchronization of qubits with the 
driven resonator takes place and the unperturbed spec- 
tral lines are replaced by one dominant spectral line at 
the driving frequency with v = uj (Fig. la,b). A similar 
phenomenon takes place for g > gdf) when the strength 
of resonator driving / oc is increased (Fig. lc,d). 

Indeed, with the growth of Up the number of photons in 
the resonator increases that leads to a stronger coupling 
between photons and qubits and eventual synchroniza- 
tion. The synchronization of qubits with the resonator 
is also clearly seen from Fig. le,f, where the spectral 
line of S{u) follows firmly the variation of driving fre- 
quency CO. The striking feature of Fig. 1 is that even 
rather different qubits with significant frequency detun- 
ings 1^^1,2 — ^o| ^ A become synchronized due to their 
coupling with the resonator getting the same lasing fre- 
quency v ^ UJ. 

For a better understanding of this phenomenon we an- 
alyze the dependence of the average number of photons 
in the resonator < n > on driving frequency uj (see Fig. 2, 
top panels). It has three pronounced maxima which up 
to quantum fluctuations correspond to three values of 
the total spin component 5'^ = Tr(p(cri"^^ + cri^^)/2) be- 
ing close to the values Sz = —1,0, 1 with the total spin 
5* = 1 (triplet state) as it is clearly seen from the data 
shown in the middle panels of Fig. 2. The whole depen- 
dence of < n > on a; is well fitted by the resonance curves 
= ripX^ /[4:{uj — uJo — AuJsJ'^-\-)^'^] where the frequency 
shift Aus^ appears due to the effective Rabi frequency 
Qr which gives oscillations between the resonant states. 
The value of Qr is induced by the coupling between 
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FIG. 2: (Color online) The average quantities of number 
of photons in the resonator < n > (top), total spin po- 
larization Sz (middle), and the concurrence of two qubits 
C (bottom) vs. the driving frequency 00. Left: identical 
qubits with Qi/cJo = ^2/(^0 = 1-2, right: different qubits 
with Qi/ujo = 1.1, ^2/(^0 = 1-2. Other parameters are 
X/ooo = 0.02, Up = 15, 5^ = 0.04. Symbols mark the values of 
the total spin Sz'- Sz > 0.2 (red/gray +), \Sz\ < 0.2 (black 
dots), Sz < —0.2 (blue/black x). Dashed curves on top panels 
show the resonance dependence ns^ (<^) (see text) with the res- 
onance shift AujSz — l-S5Szgoi>oy/ (1— < Sz >^)/(^Sz + 1). 



qubits and photons in Eq.(l) [6, 7, 8] and can be approxi- 
mated obtained as an average value of coupling that gives 
Qr = 2aSzguoy^{l- < >^){ns^ + 1) - ^S^a^Ro 
with Qro = gcooy/ns^ + 1. This gives the frequency shift 
Acos^ = dflR/dns^ which determines the resonant de- 
pendence ns^{oo) in a self consistent way. Such a theory 
gives a good description of numerical data as it is shown 
in Fig. 2 (top panel) where the corresponding values of 
< 5'^ > are taken from the middle panel. The numerical 
coefficient a smoothly varies between 1.14 and 1.57 for 
0.01 < g < 0.06. This resonance dependence is similar 
to the one qubit case discussed in [13] but the effects of 
quantum fluctuations are larger due to mutual effective 
coupling between qubits via the dissipative resonator. On 
the basis of these estimates it is natural to assume that 
the quantum synchronization of a qubit with a driving 
phase takes place under the condition that the detuning 
is smaller than the typical value of Rabi frequency 



1^1,2 -uj\< ftRo = guJo\/ns, + 1 ^ guJo\/np + 1. (3) 

This criterion assumes a semiclassical nature of photon 
field with the number of photons Up > 1. Its structure is 
similar to the classical expression for the synchronization 
tongue which is proportional to the driving amplitude 
being independent of dissipation rate [14]. The relation 
(3) determines the border for quantum synchronization 
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FIG. 3: (Color online) Time evolution of system quantities: 
average photon number < n > (top) ; components of total spin 
polarization Sx,y,z (middle) where x, 2/, 2;- components corre- 
spond to red, green, blue colors (curves from top to bottom at 
Lut/2n = 1.5x10^ (left) andc(;t/27r = 10^ (right)); concurrence 
C (bottom). Left: two identical qubits with ^1,2/000 = 1-2, 
right: two different qubits with Qi/uq = 1.1 and Qi/uo = 1.2. 
Data are shown at stroboscopic moments of time with driving 
phase (p = cot (mod 27r) = 0, Here 00/000 = 1, other parame- 
ters are as in Fig. 2. 



gc ^ 1^1,2 - ^1 /c^oV'np + 1. 

The entanglement between qubits is characterized by 
concurrence C (see e.g. the definition in [16]). Its de- 
pendence on 00 is shown in the bottom panels of Fig. 2. 
It is striking that the concurrence C can be close to 
unity not only for identical qubits but also for different 
qubits. Qualitatively, this happens due to synchroniza- 
tion of qubits induced by the resonator driving which 
makes them "quasi-identical" and allows to create the 
entangled state with Sz = 0. 

To understand the properties of the system in a bet- 
ter way we show the time evolution of its characteristics 
along a typical quantum trajectory in Fig. 3. Average 
number of photons in the resonator < n > shows tunnel- 
ing transitions between two metastable states induced by 
quantum fluctuations (top panels). There are no transi- 
tions to the third metastable state, seeing in Fig. 2 with 
three resonant curves, but we had them on longer times 
or for other realizations of quantum trajectories (see also 
Fig. 4). The life time inside each metastable state is of 
the order of thousands of driving periods and the change 
of < n > is macroscopically large (about a ten of pho- 
tons). The transition leads also to a change of total spin 
polarization components Sx,y,z = T^^{p{o-^,y,z-\-o'l^,y,z)/2) 
(middle panels). It takes place on a relatively short time 
scale ~ 1/A. The transition also generates emergence 
or death of concurrence C (or entanglement) which hap- 
pens on the same time scale 1/A (bottom panels). Nat- 
urally, C is maximal when Sx ^ 0. Remarkably, during 
long time intervals C remains to be close to its maxi- 
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FIG. 4: (Color online) Rows from top to bottom: average 
photon number (n) vs. Sz'-, concurrence C vs. Sz] resonator 
phase 00, and qubits phases 01,2 vs. the driving force phase 
(p. Columns from left to right: (i) two identical qubits with 
= 1-2, g = 0.04, (ii-iii) two different qubits, with 
Qi/cvo = 1.1, = 1.2, g = 0.04 and g = 0.001, all 

tree columns correspond to driving frequency oo/ooq = 1.01. 
The last column (iv) shows data for two different qubits with 
^i/ooo = 1.1, Q2/000 = 1.2, g = 0.04 but the driving frequency 
oo/ooo = 1 (c/. (ii)). Colors of points are explained in the text. 
Other parameters of simulations are: X/ooo = 0.02, Up = 15. 



mal value C = 1 even for the case of different qubits. 
We attribute this phenomenon to synchronization of two 
qubits by driven resonator. 

To display and characterize this phenomenon in more 
detail we determine the phases of oscillator and qubits 

via relations 00 = arctan(Im(a)/Re(a)), 0i,2 = arctan(< 

(1 2) (1 2) 

Sx > / < Sy >) respectively. The variation of these 

phases with the driving phase cp = cot (mod 27r) is shown 
in Fig. 4 for different values of system parameters corre- 
sponding to columns (i, ii, iii, iv). For identical qubits 
(column (i)) the numerical data obtained along one long 
quantum trajectory form three well defined groups of 
points in the plane of < n > and Sz (data are taken 
at stroboscopic moments of time t with a certain fre- 
quency comparable but incommensurate with oo to sweep 
all phases cp). It is convenient to mark these groups by 
three different colors corresponding to small (blue or dark 
grey), medium (black) and large (red or light grey) val- 
ues of < n > (the groups are also marked by horizon- 
tal lines). Such a classification shows that three groups 
have not only distinct values of < n > but also three 
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distinct locations in concurrence C and spin Sz- Also 
these groups show three lines in the phase plane for os- 
cillator {(po^if) and for each qubit (0i,2 5^)- Of course, 
due to quantum fluctuations there are certain fluctua- 
tions for qubit phases but the linear dependence between 
phases is seen very clearly, thus showing the quantum 
synchronization of system phases with the driving phase 
(p. Physically, the three groups correspond to the three 
triplet states of total spin S = 1. Indeed, for identical 
qubits the states with total spin values S = 1 and = 
are decoupled and the dynamics of S = state is trivial 
(see Eq. (1)). For different qubits (columns (ii), (iv)) the 
quantum synchronization between phases is also clearly 
seen even if two qubits have rather different frequency 
detunings. For the case (ii) the concurrence is smaller 
compared to the case of identical qubits (i) but by a 
change of driving frequency lu it can be increased (see the 
case (iv)) to the values as high as for identical qubits in 
(i). We note that for different qubits (e.g. for (ii)) there 
is a visible splitting of the middle group of black points 
in (< n >,Sz) plane which corresponds to states with 
mixed components of total spin S = 1 and = 0: in- 
deed, for different qubits the coupling between the states 
S = 1 and = is nonzero and such transitions can take 
place (however, on the phase planes the splitting of black 
points is too weak to be seen in presence of quantum fluc- 
tuations). Of course, the numerical data show the pres- 
ence of quantum fluctuations around straight lines in the 
phase planes. Nevertheless, this regime of quantum syn- 
chronization at g > gc is qualitatively different from the 
regime below the synchronization border g < gc where 
the points are completely scattered over the whole phase 
plane (column (iii)). In this regime (iii) the qubits rotate 
independently from the resonator which stays at flxed 
number of photons. In contrast, for g > gc two qubits 
move in quantum synchrony during a large number of 
oscillations being entangled. A single superconducting 
qubit lasing has been already achieved in experiments [9] . 
Our theoretical studies show that such two qubits, which 
in practice are always non-identical, can be made entan- 
gled and produce lasing in synchrony with each other. 
Being entangled such qubits can radiate entangled pho- 
tons in a microwave range. Therefore, the experiments 
similar to [9] but with two single-atoms lasing would be 
of great interest for entangled microwave photons gener- 
ation. 

In conclusion, our numerical simulations show that 
even two different superconducting qubits can move in 
quantum synchrony induced by coupling to a driven dis- 
sipative resonator, which can make them entangled. Such 
entangled qubits can radiate entangled microwave pho- 



tons that opens interesting opportunities for wireless en- 
tangled communication in a microwave domain. 

The work is funded by EC project EuroSQIP and RAS 
grant " Fundamental problems of nonlinear dynamics" . 



[1] A. Wallraff, D. 1. Schuster, A. Blais, L. Frunzio, 
R. S. Huang, J. Majer, S. Kumar, S. M. Girvin and 
R. J. Schoelkopf, Nature 431, 162 (2004). 

[2] J. Majer, J. M. Chow, J. M. Gambetta, J. Koch, 
B. R. Johnson, J. A. Schreier, L. Frunzio, D. I. Schus- 
ter, A. A. Houck, A. Wallraff, A. Blais, M. H. Devoret, 
S. M. Girvin and R. J. Schoelkopf, Nature 449, 443 

(2007) . 

[3] J. M. Fink, M. Goppl, M. Baur, R. Bianchetti, P. J. Leek, 
A. Blais and A. Wallraff, Nature 454, 315 (2008). 

[4] S. Filipp, P. Maurer, P. J. Leek, M. Baur, R. Bianchetti, 
J. M. Fink, M. Goppl, L. Steffen, J. M. Gambetta, 
A. Blais and A. Wallraff, arXiv:0812.2485[cond-mat] 

(2008) . 

[5] J. M. Fink, R. Bianchetti, M. Baur, M. Goeppl, L. Stef- 
fen, S. Filipp, P. J. Leek, A. Blais and A. Wallraff, 
arXiv:0812.2651 [cond-mat] (2008) . 

[6] E. T. Jaynes and F. W. Cummings, Proc. IEEE 51, 89 
(1963). 

[7] M. Tavis and F. W. Cummings, Phys. Rev. 170, 379 
(1968). 

[8] M.O. Scully and M.S. Zubairy, Quantum optics, Cam- 
bridge Univ. Press, Cambridge (1997). 

[9] O. Astafiev, K. Inomata, A. O. Niskanen, T. Yamamoto, 
Yu. A. Pashkin, Y. Nakamura and J. S. Tsai, Nature 449, 
588 (2007). 

[10] A.N. Korotkov, Phys. Rev. B 60, 5737 (1999); ibid. 67, 

235408 (2003); A.N. Korotkov and D.V. Averin, ibid. 

64, 165310 (2001). 
[11] J. Hauss, A. Fedorov, C. Hutter, A. Shnirman, and 

G. Schon, Phys. Rev. Lett. 100, 037003 (2008). 
[12] J. Gambetta, A. Blais, M. Boissonneault, A. A. Houck, 

D. I. Schuster and S. M. Girvin, Phys. Rev. A 77, 012112 

(2008). 

[13] O.V.Zhirov and D.L.Shepelyansky, Phys. Rev. Lett. 100, 
014101 (2008). 

[14] A. Pikovsky, M. Rosenblum, and J. Kurths, Synchroniza- 
tion: A Universal Concept in Nonlinear Sciences, Cam- 
bridge Univ. Press (2001). 

[15] S. Andre, V. Brosco, A. Shnirman, and G. Schon, 
arXiv:0807.4607v2 (2008). 

[16] T.E. Tessier, 1. H. Deutsch, A. Delgado, and I. Fuentes- 
Guridi, Phys. Rev. A 68, 062316 (2003). 

[17] A. Retzker, E. Solano, and B. Reznik, Phys. Rev. A 75, 
022312 (2007). 

[18] U. Weiss, Dissipative quantum mechanics. World Sci., 

Singapore (1999). 
[19] T.A. Brun, I.C. Percival, and R. Schack, J. Phys. A 29, 

2077 (1996); T.A. Brun, Am. J. Phys. 70, 719 (2002). 



